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Major advances in the observation and theory of cosmic microwave background 
anisotropies have opened up a new era in cosmology. This has encouraged the hope 
that the fundamental parameters of cosmology will be determined to high accuracy 
in the near future. However, this optimism should not obscure the ongoing need for 
theoretical developments that go beyond the highly successful but simplified stan- 
dard model. Such developments include improvements in observational modelling 
(e.g. foregrounds, non-Gaussian features), extensions and alternatives to the sim- 
plest inflationary paradigm (e.g. non-adiabatic effects, defects), and investigation of 
nonlinear effects. In addition to well known nonlinear effects such as the Rees-Sciama 
and Ostriker-Vishniac effects, further nonlinear effects have recently been identified. 
These include a Rees-Sciama-type tensor effect, time-delay effects of scalar and tensor 
t^j- ' lensing, nonlinear Thomson scattering effects and a nonlinear shear effect. Some of 

OO ■ the nonlinear effects and their potential implications are discussed. 

o 
o . 

OS . 1. INTRODUCTION 

OS 

In 1948, Alpher and Herman predicted a cosmic microwave background (CMB) radiation, 
which remained a theoretical possibility until the serendipitous discovery by Penzias and 
Wilson in 1965. A potential crisis for cosmology, arising from the apparent perfect isotropy of 
the CMB, was turned into a dramatic success when in 1992 the COBE satellite detected small 
anisotropies at the level predicted by the standard inflationary cold dark matter model. The 
confluence of theoretical and observational successes has opened up a new era of precision- 
^ ■ tested cosmology, exemplified by the upcoming launch of the MAP and Planck satellites 
(see e.g. [Q). Together with parallel advances in galactic, supernovae, lensing and other 
observations, this has promoted the hope that the fundamental parameters of the standard 
cosmological model (f2 m , h, n, . . . ) will be determined to high accuracy, and that in 
some sense, cosmology will be "solved" up to details of fine-tuning. 

However, physics rarely turns out be as simple as this, and optimism needs to be tempered 
by a realisation that future observations could entail not only the resolution of old problems, 
but also the generation of new and unexpected ones (see also 0). The new "golden age" of 
cosmology is perhaps better seen as the end of the beginning, rather than the beginning of 
the end. 
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It is clearly necessary to refine and develop the predictions of the standard model in 
readiness for the emerging observational data, and much effort has been put into this (see 
e.g. 0- [0). But it is also necessary to develop beyond the simplified standard model, in 
recognition of the rich complexity of the universe, and in anticipation of new surprises from 
future observations. This encompasses a range of theoretical advances: 

• The modelling of observations needs significant development in relation to foreground 
extraction, non-Gaussian features, data analysis, etc. (see e.g. [|]|§). 

• Extensions and alternatives to the simple inflationary paradigm need to be explored, 
including investigation of non-adiabatic effects (during inflation and reheating), of 
string cosmology (see e.g. fllPp, and of defect theories (see e.g. []TT| , p~2|| ) . Such explo- 



ration can lead to unexpected results; for example, it has recently been shown [13| 



that resonant amplification of metric perturbations in preheating after inflation could 
re-process the primordial power spectrum, and leave an imprint on CMB anistropies. 

• Nonlinear effects in CMB anisotropies require further study. 

This article focuses on the latter aspect, in particular on recently identified nonlinear 
effects. Well-known nonlinear effects include the Ostriker-Vishniac |TJJ and kinetic Sunyaev- 



Zeldovich [15] effects, which are local effects that arise from the modulation of the Doppler 
effect with inhomogeneities in the optical depth. These effects could restore small-scale 
power via hot cluster gas or re-ionisation (see e.g. [[RJ). The Rees-Sciama effect |T7| arises 



as a second-order contribution to the integrated Sachs- Wolfe term. Nonlinear integrated 
effects can in principle be significant because of the long photon paths after last scatter- 
ing. In practice the Rees-Sciama effect is typically very small. More recently, it has been 
shown [Tj| that a larger effect follows from gravitational lensing of CMB photons by density 



perturbations. 

Considerable effort has been put into calculating the imprint of these particular nonlinear 
effects on CMB anisotropies, and this work is of great importance for isolating the secondary 
anisotropies and revealing the primordial spectrum [|IIJ . What has received less attention is a 
systematic approach to nonlinear effects in general. This is important because it could reveal 
new sources of secondary anisotropies that may be significant. Furthermore, a systematic 
analysis of second-order effects forms the basis for estimating the theoretical errors associated 
with the linear analysis on which CMB science rests. 

Mollerach and Matarrese |2IJ recently developed a systematic second-order Sachs- Wolfe 



analysis, building on previous work by Pyne and Carroll They found the second- 



order metric perturbations to a matter-dominated flat universe in the Poisson gauge by 
transforming the known synchronous gauge solutions, using the methods of Bruni et al. 



22 . Then these solutions were used to find the second-order Sachs- Wolfe effect. The main 



relevant effects which they identified are: 

• in addition to the Rees-Sciama effect, another nonlinear integrated Sachs- Wolfe effect 
representing corrections to the linear gravitational wave contribution; 

• in addition to gravitational lensing by linear density perturbations, gravitational lens- 
ing by linear gravitational wave modes; 
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time delay effects of scalar and tensor lensing; 

a coupling of the velocity at last scattering with the perturbed photon wave vector. 



An alternative approach to nonlinear effects is developed in [|23|| , and will be briefly de- 
scribed here. Instead of starting from a background model and perturbing, this approach 
starts from the fully nonlinear model of the inhomogeneous universe. The approach is thus 
more general, although this generality comes at the cost of greater difficulty in quantifying 
the effects on CMB anisotropies. More importantly, the alternative approach covers non- 
linearity not only in the metric, but also in Thomson scattering, and in relative motions 
of the particle species. (Second-order scattering effects have previously been investigated 



via a different approach by Hu et al. [2~I|.) It thus offers a foundation for a comprehensive 
second- and higher-order analysis, taking into account all sources of nonlinearity. (Nonlinear 
polarisation has not been treated, although key elements of such a treatment can be found 
in Challinor's linear analysis [p5fl .) 

The key local effects identified in the new approach are: 

• at each multipole order £, a coupling of baryonic bulk velocity to the radiation bright- 
ness multipoles of order £ ± 1; 

• a coupling of the acceleration, vorticity and shear at order I to the brightness multipoles 
of order £,£±1,£±2; 

• for I 3> 1, the shear coupling is dominant and could be important. 

These local effects need to be integrated to evaluate their impact on CMB anisotropies, 
but a qualitative understanding has been achieved as a first step. 



2. NONLINEAR DYNAMICS 



The covariant Lagrangian, or 1+3 covariant, approach is based on a physical choice of 
a 4- velocity vector field u a . All variables are in principle physically measurable by co mov- 
ing observers. This approach is inherently nonlinear. It starts from the inhomogeneous 
and anisotropic universe, without a priori restrictions on the degree of inhomogeneity and 
anisotropy, and then applies the Friedmann limit when required. The basic theoretical in- 
gredients are: the covariant Lagrangian hydrodynamics of Ehlers and Ellis [261,271, and the 



perturbation theory of Hawking [2lJ and Ellis and Bruni [29] which is derived from it; the 



covariant Lagrangian approach to kinetic theory of Ellis et al. |30| ; and the covariant analysis 
of temperature anisotropies introduced by Maartens et al. [|3TJ and developed by Challinor 

331. 



and Lasenby |3^J25[[ and Gebbie et al. 

The projected symmetric tracefree (PSTF) parts of vectors and rank-2 tensors are 

V( a ) = h a b Vb, S(ab) = \h{ a c hb) d — \h cd h ab } S cd , 

where h ab = g ab + u a u b is the projector, with g ab the spacetime metric. The skew part of 
a projected rank-2 tensor is spatially dual to the projected vector S a = \e abc S^ bc \ where 
£ a bc = VabcdU d is the projection of the spacetime alternating tensor. Any projected rank-2 
tensor has the irreducible decomposition 
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Sab — ^Sh a f) + B abc S C + S( a b) , 

where S = S cd h cd is the spatial trace. A covariant vector product and its generalization to 
PSTF rank-2 tensors are 

[V, W] a = BabcV b W C , [S, Q} a = 6 a bcS b d Q Cd . 

The covariant derivative V a produces time and spatial derivatives 

j a -... b = u'T, ./"-..,, , I), ./"-..,, = //,''//", ■ • • h b f V d J e -... f . 



The projected derivative D a further splits irreducibly into a spatial divergence and curl [34 

div V = B a V a , (dwS) a = D b Sabi 

curl V a = £abcD b V c , curl S ab = e cd{a 'D c S b) d . 

Relative motion of comoving observers is encoded in the kinematic quantities: the ex- 
pansion 6 = D a u a , the 4-acceleration A a = u a = A( a ), the vorticity u a = — |curlu a , and 
the shear a ab = D( a u b y The dynamic quantities describe the sources of the gravitational 
field: the (total) energy density p = T ab u a u b , isotropic pressure p = ^h ab T ab , energy flux 
q a = —T(a)bU b , and anisotropic stress 7r ab = T< a &)> where T ab is the total energy-momentum 
tensor. The locally free gravitational field, i.e. the part of the spacetime curvature not di- 
rectly determined locally by dynamic sources, is given by the Weyl tensor C abcd . This splits 
irreducibly into the gravito-electric and gravito-magnetic fields 

E ab = C acb dU c u d = E^) , H ab = \e ac dC cd be u e = H( ab ) , 

which provide a covariant Lagrangian description of tidal forces and gravitational radiation. 

The Ricci identity for u a and the Bianchi identities V d C abcd = V[ a (— R b ] c + ^Rg b ] c ) 
produce the fundamental evolution and constraint equations governing the above covariant 
quantities [|(||7|]. Einstein's equations are incorporated via the algebraic replacement of the 



Ricci tensor R ab by T ab — \T c c g ab . These equations, in fully nonlinear form and for a general 
source of the gravitational field, are: 
Evolution: 



p + (p + p)Q + div q = -2A a q a - a ab n ab , (1) 

§e 2 + i ( 



e + ±6 2 + \{p + 3p) - div A = -a ab a ab 



+ 2uj a uj a + A a A a , (2) 
<j( a ) + f @<?a + (p + p)A a + D a p + (diwr) a = -cr ab q b 

+ [Uj,q} a - A b 7T abj (3) 

U(a) + |e^ a + |curl A a = a ab u b , (4) 

&(ab) + |©0"a6 + Eab ~ \^ab ~ D(a^b> = ~ Vc{a&b) C 

~ ^(a^b) + A {a A b ) , (5) 

E {ab) + QE ab - curl H ab + \{p + p)a ab + ±rt {ab) + |07r a6 
+ \V( a q b) = -A {a q b) + 2A c e cd(a H b) d + 3a c{a E b) c 

— ^ C £cd(aE b ) d — \(T C { a ^b)c ~ \u l '£ cd{a^ 6) , (6) 

H {ab) + QH ab + cmlE ab - ±curl7r ab = 3a c{a H b) c - u c e cd{a H b) d 

— 2A c s cd i ya E b d — ^uJ( a q b ) + \ct' \ a e b ) cd q d . (7) 



Constraint: 



divto = A a u a , (8) 

(di V( x) a - curler - |D a 9 + g a = -2[u,A] a , (9) 

curler^ + D {a u b) - H ab = -2A {a u b) , (10) 
(divE) a + |(div7r) a - §D a p + §6g a = [a,H] a 

-3H ab u b + ±a ab q b -l[uj,q] a , (11) 
(divH) a + |curlg a - (p + p)u a = -[a,E) a 

-\[a^} a + 3E ab u b -\u ab uj b . (12) 



If the universe is almost Friedmann, then quantities that vanish in the Friedmann limit 
are 0(e), where e is a dimensionless smallness parameter, and the quantities are suitably 
normalised (e.g. \J a ab a ab / Q < e, etc.). Linearisation reduces all the right hand sides of the 
evolution and constraint equations to zero. 

For a given choice of fundamental frame u a , each of the species I which source the 
gravitational field has relative velocity vj and 4-velocity 

v? = >y I (u a + v?) , v?u a = 0. (13) 

If the 4-velocities are close, i.e. if the frames are in non-relativistic relative motion, then 
0(v 2 ) terms may be dropped from the equations, except if we include nonlinear kinematic, 
dynamic and gravito-electric/magnetic effects, in which case, for consistency, we must retain 
0(e°v 2 ) terms such as pv 2 , which are of the same order of magnitude in general as 0(e 2 ) 
terms. This is the physically relevant nonlinear regime, i.e. the case where only nonrelativis- 
tic bulk velocities are considered, but no restrictions are imposed on non-velocity terms, and 
we neglect only terms 0(ev 2 ,v 3 ). Second-order effects involve neglecting also terms 0(e 3 ), 
while linearisation means neglecting terms 0(e 2 ,ev,v 2 ). 

The dynamic quantities in the evolution and constraint equations ([ID— ( |P2l) are the total 
quantities, with contributions from all dynamically significant particle species. Thus 

T ab = T ab = pu a u b + ph ab + 2q {a u b) + 7i ab , (14) 
i 

T ab = Pl uX + hhf + 2q$ a vP + < 6 • (15) 

The dynamic quantities p n ■ ■ ■ in equation (|l^) are as measured in the /-frame. For cold 
dark matter (CDM), baryons, photons and neutrinos 

Pa =0 = <g = n« b , g ; = = <\ p R = \p R , Pn = 1 Pn , 

where we have chosen the unique 4-velocity in the CDM and baryonic cases which follows 
from modelling these fluids as perfect. The cosmological constant is characterized by 

P V = -P V = -A, < = = <\ < = 0. 

The conservation equations for the species are best given in the overall ii a -frame, as are 
the evolution and constraint equations above. This requires the expressions for the partial 
dynamic quantities as measured in the overall frame. The inverse velocity relation is 
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where v* a u ia = 0, and v* a v* a = vfv ia . Then the dynamic quantities of species I as measured 
in the overall u a -frame are 



p; 

q *a 



71 



*ab 



ft + {iftf (ft + ft) + 2 7/ g>,a + <Vft&} , 
^ + I bS (Pi + ft) + H<l>ia + <Vft&} > 

?;+(ft +ftK 

+ {{ir - 1)?; - 7^V^ a + 7> 7 2 (a + Pj) < 

+ TifV - 7T* c v lb v lC u a ] , 
^ + {-2u^ c v lC + ir; c v lb v lC u a u b } 
+ {-l^c^d/i 06 + 7/ (ft + ft) ^ (a ft 6) + 27,^^} 



(16) 
(17) 



(18) 



(19) 



These are the nonlinear generalisations of well-known linearised results, which correspond to 
removing all terms in braces, dramatically simplifying the expressions. To linear order, there 
is no difference in the dynamic quantities when measured in the /-frame or the fundamental 
frame, apart from a simple velocity correction to the energy flux. But in the general nonlinear 
case, this is no longer true. 

The total dynamic quantities are simply given by 



P 



Eft*, p = £ft% <f = EC 



71 



*ab 



A convenient choice for each partial 4-velocity u a is the energy frame, i.e. q a = for each 
/. In the fundamental frame, the partial energy fluxes do not vanish, i.e. q* a 7^ 0, and the 
total energy flux is given by 

<f = E [(ft + ft) K + <V + °K> ft 3 / • 

Then the dynamic quantities of CDM as measured in the fundamental frame are 



* 2 * 1 2 2 
P c =7 c ft, P = H% v cPci 



^=%Pc%, 



71 



while for baryonic matter: 



D 



ll (l + w B vl) ft , p* B 



* ab -=^0 v {a v b) 



7 S 2 (1 + w B )p B v a B , t£* = 7fl 2 (l + w B )p B v^ , 



(20) 
(21) 



(22) 
(23) 



where w B = p B /p B - For radiation and neutrinos, the dynamic quantities relative to the 
ii a -frame are found directly via kinetic theory below. 

The total energy-momentum tensor is conserved, i.e. VbT ab = 0, which is equivalent to 
the evolution equations ([!]) and (|3|). The partial energy-momentum tensors obey 



v T ab 



J* = U*u a + M* a , 



(24) 
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where U* is the rate of energy density transfer to species / as measured in the n a -frame, 
and M* a = M*( a < is the rate of momentum density transfer. Thus 

j c a = o = j;, r R =-r B =u T u a + M^ 

where the Thomson rates are, to 0(ev 2 ,v B ), 

K = w (IpX - € + *, 



as shown below. Here n E is the free electron number density, and a T is the Thomson cross- 
section. Note that beyond linear order there is energy transfer, i.e. U T 7^ 0. 
Using equations ©-(|23]) in @§, we find that, to 0(ev 2 ,v 3 ), for CDM 



Pc + Q Pc + Pc div v c = ~ (Pc v c) ~ i v c Q Pc 

-v«D aPc -2p A a v«, (25) 

t + l Qv c + A<1 = - ff >c + k VcY - v b c V h v a c , (26) 

and for baryonic matter 

P B + ©(! + w b)Pb + (1 + w B )p B divv B = - [(1 + w B )p B v 2 

- §t£e(i + w b ) Pb - <d„ [(1 + w b ) Pb ] 

- 2(1 + w B )p B A a v a B - n E a T (f p> 2 - q* R a v Ba ) , (27) 
(1 + w B )v a B + (§ - c 2 ) 9t£ + (1 + w B )A a + p B l Wp B 

+ Pb^ e °t {pX - c) = (1 + w B )A bV y 

- (1 + w B )a\v b B + (1 + w B )[u, v B f - (1 + w B )v b D b v a B 

+ 4(1 + w B )(divv B X ~ P B l n E a T ^\ b , (28) 

where c 2 = p B /p B (this equals the adiabatic sound speed only to linear order). These are the 
nonlinear energy conservation and relative velocity equations for matter. Linearization re- 
duces the right hand sides to zero, dramatically simplifying the equations. The conservation 
equations for the massless species (radiation and neutrinos) are found below. 

3. NONLINEAR THOMSON SCATTERING 

In covariant Lagrangian kinetic theory |30| the photon 4-momentum is split as 



p a = E{u a + e a ) , e a e a = 1 , e a u a = 



where E = —u a p a is the energy and e a = p^ / E is the direction, as measured by a co- 
moving (fundamental) observer. The photon distribution function is expanded in covariant 
harmonics 
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f{x,p) = f{x, E,e)=F + F a e a + F ab e a e b + 



J2F Al (x,E)e^\ 



where e Ae = e ai e a2 ■ ■ ■ e ae and 



Fa-b — F( a ... b ) F a ... b — F( ffl ...6) , F a ... b u b — — F a ... bc h bc . 



The first 3 multipoles arise from the radiation energy-momentum tensor, 

B {x)=JPP f{x,p)d p = p R u u + gp^/i + 2q R y u' + n R 
The radiation brightness multipoles are 



n 



E S F dE 



so that (dropping asterisks) II = p R /4ir, U a = 3q R /47i and n ab = 157r^ fc /87r. These multipoles 
define the temperature fluctuations [pT|] . 
The Boltzmann equation is 



df 



df 



dv ^ dx a 



" ' T a bc p b p c0j 



C[f] = b + b a e a + b ab e a e b + 



dp 



where the collision term C[f] determines the rate of change of / due to emission, absorption 
and scattering processes. For simplicity, the effects of polarization are neglected (see for 
a linear treatment), and 

C[f] = a T n E E B [f(x,p) - f(x,p)] , 
where E B = —p a u R is the photon energy relative to the baryonic frame u a , and |32 

Jf(x,p') [l + (eXa) 



f(x,p) 



16tt 



dn' B . 



(29) 



Here e' „ is the initial and e a is the final direction, so that 



P 



P" 



E B [u a + e a 



where we have used E' — E B , which follows since the scattering is elastic. The exact forms 
of the photon energy and direction in the baryonic frame are 



E% l~v a e a 



e a + ll(vie b -vl)u a + 1 l(v b e b -l)v: 



E n 



B 7 s (l-^e c ) 

Since the baryonic frame will move non-relativistically relative to the fundamental frame 
in all cases of physical interest, it is sufficient to linearize only in v B , and not in the other 
quantities. Thus we drop terms in 0(ev 2 , v 3 ) but do not neglect terms that are 0(e°v 2 , ev B ) 
or 0(e 2 ) relative to the background. It follows from equation (p9|) that 
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4tt / fE*dE B = {p R ) B + f « ) B e Ba e Bb 



where the dynamic radiation quantities are evaluated in the baryonic frame. Transforming 



back to the fundamental frame, we find, to 0(ev E ,v E ), 



{P R ) B = P R I 1 + Kj - 2 C*V, 
(-n ab ) = 7c ab + 2v RC 7T c{a u b) - 2q {a v b) + ±p R v {a v b) . 

^ R >B R BC r ^R B 3 r R B B 



Using the above equations and various identities p3|, and defining the energy- integrated 
scattering multipoles 



K Ae = J E 2 b Ae dE. 



we find that to 0(ev E ,v E ) 



K = n E a T 
K a = -n E a T 



K 



ab 



K 



abc 



and, for £ > 3: 



-T\v - -Wv 

3 iiy fl 3 ii u Ba 



IF - AIiv a - f n a V 



9 TTafc 
10 



n 



abc 



jn^) - lu abc \ d 



3Uv {a v b) 



(30) 
(31) 
(32) 
(33) 



-n E a T 



6 + 1 1 u Ata v 
21 + 3 Ba 



(34) 



Equations (|3~0"D -( j34]) are a nonlinear generalisation of the linearised Thomson scattering re- 
sults [j32|| . They show clearly the coupling of baryonic bulk velocity to the radiation multipoles, 
arising from local nonlinear effects in Thomson scattering. 

The multipoles of E~ 1 df/dv are derived in |30[ and f23|, using different methods. The 
result is 



F {Al) -\QEF' Ai + V {ae F Al _ l) + 



(*+l) 
(2£ + 3) 



D a F, 



(2£ + 3) 
2^ 



E v ~ l F, 



{£ + !){£ + 2) 
(2£ + 3)(2£ + 5) 



171-1/2 7713/2 771 

(27T^) E r F ^ 



{At 



E t+3p 



A 



at) 



abAf 



a 



E 2 -«E 



a e _ 1 a t ) j 



where a prime denotes d/dE. This result is exact and holds for any photon or (massless) 
neutrino distribution in any spacetime. Integrating, it leads to 
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K Af = u {Ae) + ien Af + D (af n 



(a e 11 -A e _ 1 } 



+ 



+ D 



■D b n 



bA, 



A b Tl, 



bA, 



(2£ + 3) 
(£-!)(£ +!)(£ + 2) b 

(2£ + 3)(2£ + 5) bc ^ + 

(£ + 2)a (a , a ,_ i n^_ 2> . 



(2£ + 3)" 

(£ + 3)A {ae U Ai _ l} - ^A^IU^/ 
U 



{21 + 3) 



a 



i>6 



(35) 



For decoupled neutrinos, K At = 0, while for (unpolarised) photons undergoing Thomson 
scattering, the left hand side of equation ( p5|) is given by equations (|3f]D-([3"lD , which are exact 
in the kinematic and dynamic quantities, but first order in the baryonic bulk velocity. These 
equations show the coupling of acceleration, vorticity and shear to the radiation multipoles 
that arises at nonlinear level. 

The monopole and dipole of equation (|35| ) imply photon conservation of energy and 
momentum density (to 0(ev\vt)): 



Pr + |©P« + + 2A °€ + = U e°t (lP R vl - <t> B a) , 



g<°> + f 0<£ + lp R A a + \V a p R + D< + a\q b R - [ U , q R ] a + A b n 



ab 
R 



n, 



a f i a _ a , ab 



%b 



(36) 
(37) 



The quadrupole evolution equation is 



97r a6 + ^p R a afc + |D<Y> + — D c n 

R 15 r R 5 fl 35 



rabc 



(<V fe > c 
J? 



32VT -ppahcrf 

315 



(38) 



and the higher multipoles (I > 3) evolve according to 

n<^> + |en^ + d<°ti a ^ + ;£ +1 ^ D 6 n 6i 



(i + !)(*- 2) 
(2£ + 3) 



(2£ + 3)~ 

A b tt bAe + {£ + 3)A {ae U Ae ~ l} - £u b e bc {a m A ^ 



(£-l)(£ + l)(£ + 2) bcA 

a bc LL + 



(2£ + 3)(2£ + 5) 
- (£ + 2)a {aeae ~ 1 U Ae - 2) 



(2£ + 3) 



a b {a m Ae - i)b 



—n B a T 



]J A e — ]j( A i-i v a i) _ 



+ 1 



, U Aea v 

21 + 3 sa 



(39) 



For £ = 3, the second term in square brackets on the right of equation ( |39"D must be multiplied 
by §. The temperature fluctuation multipoles are determined by the radiation brightness 



multipoles Ii A [EH 



These equations show in a transparent and covariant form precisely which physical effects 
are directly responsible for the evolution of CMB anisotropies in an inhomogeneous universe. 
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They apply at second and higher-order, and are readily specialised to the linear case. They 
show how the matter generates anisotropies: directly through interaction with the radia- 
tion, as encoded in the Thomson scattering terms on the right of equations (|36|)-(|39|); and 
indirectly through the generation of inhomogeneities in the gravitational field via the field 
equations (]!]) (p~2|) and the evolution equation fl2"5|) for the baryonic velocity v®. This in turn 
feeds back into the multipole equations via the kinematic quantities, the baryonic velocity 
v®, and the spatial gradient D a p fl in the dipole equation (|57|). The coupling of the mul- 
tipole equations themselves provides an up and down cascade of effects, shown in general 
by equation (PP|). Power is transmitted to the I- multipole by lower multipoles through the 
dominant (linear) distortion term D^IT 4 *- 1 ^, as well as through nonlinear terms coupled 
to the 4-acceleration (A^U^- 1 ^), baryonic velocity (^IT 4 *-^), and shear (er< a <°'-in A *- 3 >). 
Simultaneously, power cascades down from higher multipoles through the linear divergence 
term (divIT) A *, and the nonlinear terms coupled to A a , v" and a ab . The vorticity coupling 
does not transmit across multipole levels. 

4. CONCLUSION 

A range of nonlinear effects is identified via a systematic covariant analysis. These include 
the following. 

• Nonlinear relative velocity corrections, as exemplified in the dynamic quantities in 
equations (|T6|) — (|19D and in the bulk velocity equations fl26|) and 

• Nonlinear Thomson scattering corrections affect the baryonic and radiation conserva- 
tion equations, entailing a coupling of the baryonic bulk velocity v a to the radiation 
energy density, momentum density and anisotropic stress. The evolution of the radia- 
tion quadrupole 7i^ b also acquires nonlinear Thomson corrections, which couple to 
the radiation dipole q° and octopole U abc . Linearization, by removing these terms, has 
the effect of removing the nonlinear contribution of the radiation multipoles IP 4 *'* 1 to 
the collision multipole K l . 

• Nonlinear kinematic corrections introduce additional acceleration and shear terms. 
Vorticity corrections are purely nonlinear, i.e. a linear approach could give the false im- 
pression that vorticity has no direct effect at all on the evolution of CMB anisotropies. 
However, for very high £, i.e. on very small angular scales, the nonlinear vorticity 
term could in principle be non-negligible. The general evolution equation fl39|) for the 
radiation brightness multipoles H A * shows that five successive multipoles, i.e. for £ — 2, 
• • •, £ + 2, are linked together in the nonlinear case. The 4-acceleration A a couples 
to the £ ± 1 multipoles, the vorticity u a couples to the £ multipole, and the shear a ab 
couples to the £ ± 2 and £ multipoles. All of these couplings are nonlinear, except 
for £ = 1 in the case of A a , and £ = 2 in the case of o~ a b- These latter couplings 
that survive linearisation are shown in the dipole equation (|37| ) (i.e. p R A a ) and the 
quadrupole equation (|38D (i.e. p R a ab ). The latter term drives Silk damping during 
the decoupling process. The disappearance of most of the kinematic terms upon lin- 
earisation is further reflected in the fact that the linearised equations link only three 
successive moments, i.e. £, £±1. 
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A crucial feature of the nonlinear kinematic terms is that some of them scale like £ for 
large £, as already noted in the case of vorticity. There are no purely linear terms with 
this property, which has an important consequence, i.e. that for very high £ multipoles 
(corresponding to very small angular scales in CMB observations), certain nonlinear 
terms may reach the same order of magnitude as the linear contributions. (Note that 
the same effect applies to the neutrino background.) The relevant nonlinear terms in 
Eq. @ are (for £ > 1): 

- £ [\a bc Ii bcAl + a W^-in^- 2) 
+ ^(^n Af - l} + \A b Tl bAt + Lu b e bc {ae Tl Ae - l)c 

Any observable imprint of this effect will be made after last scattering. In the free- 
streaming era, it is reasonable to neglect the vorticity relative to the shear. We can 
remove the acceleration term by choosing u a as the dynamically dominant cold dark 
matter frame (i.e. choosing = 0). It follows that the nonlinear correction to the 
rate of change of the linearised fluctuation multipoles is 

5{tl Ae ) ~ £ [\a bc Il hcAl + <7< a ^-in A *- 2) ] for £ > 1 . (40) 

The linear solutions for rU £ and a ab can be used in equation (|40| ) to estimate the 
correction to second order. Its effect on observed anisotropies will be estimated by 
integrating 5(JX Al ) from last scattering to now. 



Further quantitative analysis of the new nonlinear effects in is needed to iden- 

tify more clearly which effects could be observationally significant. This would also clarify 
the relationship between various nonlinear effects that have been derived under different 
assumptions and using different approaches [p^,p!5|,p^,p0|,p^,p4]1 . 
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